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We study the dynamics of entanglement in a two-qubit system interacting with a squeezed thermal 
bath via a dissipative system-reservoir interaction with the system and reservoir assumed to be in a 
separable initial state. The resulting entanglement is studied by making use of a recently introduced 
measure of mixed state entanglement via a probability density function which gives a statistical and 
■ geometrical characterization of entanglement by exploring the entanglement content in the various 

subspaces spanning the two-qubit Hilbert space. We also make an application of the two-qubit 
f^*) dissipative dynamics to a simplified model of quantum repeaters. 
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I. INTRODUCTION 

Oh! Open quantum systems take into account the effect of the environment (reservoir or bath) on the dynamical evolution 
■ of the system of interest thereby providing a natural route for discussing damping and dephasing. One of the first 
§ | testing grounds for open system ideas was in quantum optics [l[ . Its application to other areas gained momentum from 
the works of Caldeira and Leggett Q, and Zurek among others. The total Hamiltonian is H = Hs + Hr + Hsr, 
q-i where S stands for the system, R for the reservoir and SR for the system-reservoir interaction. Depending upon 
the system-reservoir (S — R) interaction, open systems can be broadly classified into two categories, viz., quantum 
non-demolition (QND) or dissipative. A particular type of quantum nondemolition (QND) S — R interaction is given 
by a class of energy-preserving measurements in which dephasing occurs without damping the system, i.e., where 
[Hs,Hsr] = while the dissipative systems correspond to the case where [Hs,Hsr] ^ resulting in decoherence 
along with dissipation Q . 

A prototype of dissipative open quantum systems, having many applications, is the quantum Brownian motion 
of harmonic oscillators. This model was studied by Caldeira and Leggett Q for the case where the system and its 
environment were initially separable. The above treatment of the quantum Brownian motion was generalized to the 
jhysically reasonable initial condition of a mixed state of the system and its environment by Hakim and Ambcgaokar 
s * . H|, Smith and Caldeira [||, Grabert, Schramm and Ingold 0], and for the case of a system in a Stern-Gerlach potential 
gj, and also for the quantum Brownian motion with nonlinear system-environment couplings among others. 

The interest in the relevance of open system ideas to quantum information has increased in recent times because 
of the impressive progress made on the experimental front in the manipulation of quantum states of matter towards 
quantum information processing and quantum communication. Myatt et al. [lfjj ] and Turchette et al. [TTI | have 



performed a series of experiments in which they induced decoherence and decay by coupling the atom (their system- 
s') to engineered reservoirs, in which the coupling to, and the state of, the environment are controllable. 

Quantum entanglement is the inherent property of a system to exhibit correlations, the physical basis being the 
non-local nature of quantum mechanics [12j, and hence is a property that is exclusively quantum in nature. Entan- 
glement plays a central role in quantum information theory [l3| as in interesting non-classical applications such as 
quantum computation [bij and quantum error correction [lq . A number of methods have been proposed for creating 
entanglement involving trapped atoms [l?], HH ■ 

An important issue is to study how quantum entanglement is affected by noise, which can be thought of as a 
manifestation of an open system effect (l9| . In [2p| entanglement of a two-mode squeezed state in a phase-sensitive 
Gaussian environment was studied and the criteria for the necessary and sufficient condition for separability of 
Gaussian continuous- variable states (2l| was employed as a measure of entanglement. In [22j the entanglement between 
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charge qubits induced by a common dissipative environment was analyzed using concurrence as the measure. Some 
recent experimental investigations on the influence of decoherence on the dynamics of entanglement have been made 
in [23|, [24| . In a related work [2f| , this issue was taken up with the noise coming from the effect of the environment 
modelled by a QND S — R interaction. Here we complement this program by studying the effect of noise, modelled 
by a dissipative S — R interaction with the reservoir in an initial squeezed-thermal state [U, [H| , on the entanglement 
evolution between two spatially separated (and initially uncorrelated) qubits, brought out by interaction with the 
bath. This would be of relevance to evaluate the performance of two-qubit gates in practical quantum information 
processing systems. 

Since we are dealing here with a two qubit system which very rapidly evolves into a mixed state, a study of 
entanglement would necessarily involve a measure of entanglement for mixed states. Entanglement of a bipartite 
system in a pure state is unambigious and well defined. However, mixed state entanglement (MSE) is not so well 
defined. Thus, although a number of criteria such as entanglement of formation (27l. l28l. [29| and separability [30l | exist, 
there is a realization 27[ that a single quantity is inadequate to describe MSE. This was the principal motivation for 
the development of a new prescription of MSE (3l| in which it is characterized not as a function, but as a probability 
density function (PDF). The known prescriptions such as concurrence and negativity emerge as particular parameters 
that characterize the probability density. We will principally make use of this measure in our study of entanglement 
in the two-qubit system. 

The plan of the pap er is as follows. In Section II, we recapitulate for consistency, the recently developed entanglement 
measure of MSE [3l|. In Section III, the master equation describing the dynamical evolution of the two-qubit system 
interacting with a squeezed thermal bath, via a dissipative S — R interaction, is given which is then used in Section 
IV, to study in detail the dynamics of the system interacting with a vacuum bath with zero bath squeezing in Section 
IV(A) and with a general squeezed thermal bath in Section IV(B). Section V deals with the entanglement analysis of 
the two-qubit open system using the PDF as a measure of entanglement. We compare it with the usual measure of 
MSE, concurrence. We dwell on the scenarios where the two qubits effectively interact via localized S — R interactions, 
called the independent decoherence model, as also when they interact collectively with the bath, called the collective 
decoherence model. In Section VI, we make a brief application of the model to practical quantum communication, in 
the form of a quantum repeater [32|, HH. In Section VII, we make our conclusions. 



II. CHARACTERIZATION OF MIXED STATE ENTANGLEMENT THROUGH A PROBABILITY 

DENSITY FUNCTION 

Here we briefly recapitulate the characterization of mixed state entanglement (MSE) through a PDF as developed 
in (3lj . As pointed out in the Introduction, the above criterion was evolved from the motivation that for the charac- 
terization of MSE, a single parameter is inadequate. The basic idea is to express the PDF of entanglement of a given 
system density matrix (in this case, a two-qubit) in terms of a weighted sum over the PDF's of projection operators 
spanning the full Hilbert space of the system density matrix. The PDF of a system in a state which is a projection 
operator p = jj^-m of rank M is defined as: 

jdHn M 5{£^ -£) 

Pn M {t) = 77— , (1) 

J d,Hn M 

where J dTi]j A[ is the volume measure for Wn M , which is the subspace spanned by 11m. The volume measure is 
determined by the invariant Haar measure associated with the group of automorphisms of J dHn M , modulo the 
stabilizer group of the reference state generating TLn M ■ Thus for a one dimensional projection operator, representing 
a pure state, the group of automorphisms consists of only the identity clement and the PDF is simply given by the 
Dirac delta. Indeed, if p = Hi = \ip)(ip\, the PDF has the form V p {£) = 5(£ — £^) thereby resulting in the description 
of pure state entanglement, as expected, by a single number. The entanglement density of a system in a general mixed 
state p is given by resolving it in terms of nested projection operators with appropriate weights as 

p = (Ai - A 2 )ni + (A 2 - A 3 )n 2 + (Aat_i - Ajv)nAr_i + \ N IL N 

N 

= ]T A A/ n M , (2) 

M=l 

where the projections arc Um = Y^jLi iV'jK^'jIj wrtn M = 1,2,..., A and the eigenvalues Ai > A 2 > i.e., the 
eigenvalues are arranged in a non-increasing fashion. Thus the PDF for the entanglement of p is given by 

N 

V p {£)= "MVn M (£), (3) 

M=l 
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where the weights of the respective projections Vn M (£) are given by u>m = Am/Ai. For a two qubit system, the density 
matrix would be represented as a nested sum over four projection operators, n 1; n 2 , IT3, II4 corresponding to one, 
two, three and four dimensional projections, respectively, with III corresponding to a pure state and II4 corresponding 
to a a uniformly mixed state, is a multiple of the identity operator. The most interesting structure is present in II2, 
the two-dimensional projection, which is characterized by three parameters, viz. £ CUS p, the entanglement at which 
the PDF diverges, £ m ax, the maximum entanglement allowed and ?2(£moi), the PDF corresponding to £ m ax- The 
three dimensional projection II3 is characterized by the parameter £±, which parametrizes a discontunity in the 
entanglement density function curve. By virtue of the convexity of the sum over the nested projections it can 
be seen that the concurrence of any state p is given by the inequality C p < (Ai — A 2 )Cn! + (A 2 — A 3 )Cn 2 - Thus 
while the concurrence for a three and four dimensional projection is identically zero, through the PDF one is able 
to make a statement about the entanglement content of these spaces. The fact that the PDF ^ enables us to 
study entanglement of a physical state by exploiting the richness inherent in the subspaces spanned by the system 
Hilbert space makes it an attractive statistical and geometric characterization of entanglement, of which an explicit 
illustration is made in Section V. 



III. TWO-QUBIT DISSIPATIVE INTERACTION WITH A SQUEEZED THERMAL BATH 

We consider the Hamiltonian, describing the dissipative interaction of N qubits (two-level atomic system) with the 
bath (modelled as a 3-D electromagnetic field (EMF)) via the dipole interaction as [34| 

H = Hs + Hr + Hsr 

N N 
= ^n S n + E fc^L 6 *. + 1/2) ~ ^EE^'i^X^ + S~)b n3 - h.C.]. (4) 

n=1 ks ks 11=1 

Here p n are the transition dipole moments, dependent on the different atomic positions r n and 

St = \e n )(g n \, S~ = \g„)(e n \, (5) 
are the dipole raising and lowering operators satisfying the usual commutation relations and 

Sn = \(\ e n)(e n \ - \g n )(g n \), (6) 
is the energy operator of the nth atom, while bt bt are the creation and annihilation operators of the field mode 

ks ks 

(bath) ks with the wave vector k, frequency lju and polarization index s = 1,2 with the system-reservoir (S-R) 
coupling constant 

hsifn) = {^f'^ k /^. (7) 

Here V is the normalization volume and e^ s is the unit polarization vector of the field. It can be seen from Eq. 
([7]) that the S-R coupling constant is dependent on the atomic position r„. This leads to a number of interesting 
dynamical aspects, as seen below. From now we will concentrate on the case of two qubits. Assuming separable initial 
conditions, and taking a trace over the bath the reduced density matrix of the qubit system in the interaction picture 
and in the usual Born-Markov, rotating wave approximation (RWA) is obtained as [34[ 



^ = -$r[H§,p] ~ \ E r ^ t 1 + X\(P S tSj + SfSjp - 2Sj P Sf) 

2 2 

- - J2 r i3 N( P Sr S + + SiS+p - 2S+pSH + - £ T. lJ M(pS+S+ + S+S+p - 2S+pS+) 

2 

+ - ]T Y^ipSrSj + SrS-p - 2SjpSr). (8) 



In Eq. 



N = AT th (cosh 2 (r) + sinh 2 (r)) + sinh 2 (r), (9) 
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M 



--sinh(2r)e t *(2iV th 



(10) 



with 



LO 



(11) 



and 



N th = 



e k sT - 1 



(12) 



Here iVth is the Planck distribution giving the number of thermal photons at the frequency lu and r, $ are squeezing 
parameters. The analogous case of a thermal bath without squeezing can be obtained from the above expressions 
by setting these squeezing parameters to zero, while setting the temperature (T) to zero one recovers the case of the 
vacuum bath. Eq. ([5]), for a single qubit case, can be solved using the Bloch vector formalism (cf. [35[) and also 
in the superoperator formalism [3r3 | . Here the assumption of perfect matching of the squeezed modes to the modes 
of the EMF is made along with, the squeezing bandwidth being much larger than the atomic linewidths. Also, the 
squeezing carrier frequency is taken to be tuned in resonance with the atomic frequencies. 
In Eq. ©, 



2 

71=1 



(13) 



where 



-[l-(A.r. y -) 2lC ° S(fc0 ^ ; 



k n 



+ [l-3(£.f y ) a 



sin(fc r lJ ) cos(fc ry). 



+ 



(14) 



Here p, = fii = \±2 and f%j are unit vectors along the atomic transition dipole moments and fi 



respectively. 



Also fco = ijJq/c, with ujq being as in Eq. (fTI]) . 



The wavevector ko = 27t/Aq, Aq being the resonant 



wavelength, occuring in the term fcor^ sets up a length scale into the problem depending upon the ratio r^/Ao- 
This is thus the ratio between the interatomic distance and the resonant wavelength, allowing for a discussion of the 
dynamics in two regimes: (A), independent decoherence where k^.r^ ~ > 1 and (B). collective decoherence where 
ko-Tij ~ — > 0. The case (B) of collective decoherence would arise when the qubits are close enough for them 
to feel the bath collectively or when the bath has a long correlation length (set by the resonant wavelength Ao) in 
comparison to the intcrqubit separation r,j . Qjj (|14j) is a collective coherent effect due to the multi-qubit interaction 
and is mediated via the bath through the terms 



Sirehc 3 



(15) 



The term Tj is present even in the case of single-qubit dissipative system bath interaction [35l . l36j | and is the spontaneous 
emission rate, while 



where i ^ j with 



F{k n 



Tij = Tji = x 1,1, l-ikui;,) 
sin(fcor 



konj 

cos(fc r i j) _ sin(fc r i j) 

(konj) 2 (konj) 3 



(16) 



[1 - (V- r il) l—TZ + [1 - KH-Tij) J 



(17) 



Tij (fT5j) is the collective incoherent effect due to the dissipative multi-qubit interaction with the bath. For the case 
of identical qubits, as considered here, f2i2 = ^21, Ti2 = I^i and Ti = T2 = T. 
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IV. DYNAMICS OF THE TWO-QUBIT DISSIPATIVE INTERACTION WITH A VACUUM AND 

SQUEEZED THERMAL BATH 

Here we present the solutions of the density matrix equation (J8J) for the case of a two-qubit system interacting with 
a (A), vacuum bath and (B). squeezed thermal bath. These results will be of use in the investigation of the dynamics 
of entanglement subsequently. 



A. Vacuum bath 



Here we proceed as in [34| and obtain the reduced density matrix from Eq. ([8]), setting T and bath squeezing to 
zero, and by going over to the dressed state basis, of the collective two-qubit dynamics, obtained from the Hamiltonian 



Hg ([131) as 



Iff) 
\s) 

\a) 
\e) 



m)m), 

-^(|ei)|5a) + |fli)|e 2 )), 

^(|ei)|0»>-|si>|e 2 », 
|ei}|e 2 ), 



(18) 



with the corresponding eigenvalues being E g = —Hujo, E s = Ml\2 1 E a = — ?iSli 2 and E g = Hluq. The reduced density 
matrix in the dressed state basis (fl8|) can be obtained from Eq. ([8]) as 



where 



dp i [TJ , , ,dp dp 



h[u (\e)(e\-\g}(g\) + n 12 (\s}(s\-\a}(a\ 



(19) 
(20) 



and 



[ dt h 



dp 

[ dt )a 



- 5 ( r + r i 2 ) (e| + |*>W)P + P(|e)(e| + \s)(s\) 
2(\s)(e\ + \g)(s\)p(\e)(s\ + \s)(g\)}, 



--(r - T 12 ) [(|e)(e| + |a)(o|)p + p(|e)(e| + \a)(a\) 
2(| a )(e|-| ff )(o|)p(|e)(a|-|a)(ff|)]. 



(21) 



(22) 



From the dressed state basis (fT5|) . it can be seen that the two-qubit problem can be thought of as an equivalent four- 
level system. For the case where r^/Ao — > 0, i.e., when the interatomic separation is much smaller than the resonant 
wavelength, constituting the Dicke model [33], (dp/dt) a = and the problem reduces to an effective three- level system. 
The Eq. (fl9|) can be solved to yield the various density matrix elements as follows: 



Pss(t) 



,-(r+r 12 )t 



Pss{0) 



Pee (t) = e - 2I Ve(0) 

(r + r 12 ) 



(r-r 12 ) 



(1 _ e -(r-r 12)t)e - ( r + r 12 ) Vee(0); 



(23) 
(24) 



Paa (t) = e-^-^WW + £-^4(1 - e - (r+ri2)t )e-( r - r -)Vee(0), 

(1+1 12) 



(25) 



PuS) = P S9 (0) + (l-e-( r+ri2 )*)p ss (0) + (l-e-( r " r -' t K a (0) 

"(r + r 12 ) 



2r 



i - 



1 r (r + Fl2) (i - e -(r-r 12 )t) + (£_£lHZ le -( 



l (r-r 12 ) L 2 

(r-Tiz) r , _ (r -r 12 )tN _ ( r - r i z; 
(r + r 12 )l l ' 2T 



-(i- e - 2rt )} 



2 

Pee(0). 



12)t| 



(26) 
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The Eqs. (|23|) to (|26p give the dynamics of the population of the two-qubit system interacting with a vacuum bath 
and p ee (t) + p ss (t) + p aa {t) + p gg {t) = p ee (0) + p ss (0) + p aa (0) + P gg (0)- The off- diagonal terms of the density matrix 



are: 



Pse (t) = P* es (t), (27) 
p eg (t) = e-^ ot e- rt p eg (0), 

P 9 e(t) = P* eg (t), (28) 
Pea (t) = e^'+^e-^-^peaiO), 

Pae{t) = Pl a {t), (29) 

Psait) = e- j2n - t e- rt p S a(0), 

p as (t) = p* sa (t), (30) 
n (+\ - p -i(,u) -n 12 )t -i(r-r 12 )t ^ (r-T^) i(lJJ „-n 12 )t 

Pag{t) - e e p ag (U) ( T2+m 2 2 f 

x e -^ (r_ri2) * [2fti 2 e- rt sin(2f7i 2 t) + T(l - e~ rt cos(20i 2 t))] p ea (0) 

+ l -^^-e-^-^ t e-^ r - r ^ t [2f7 12 (l- e - rt cos(2£W)) 

- Te- rt sm(2n 12 t)]p ea {0), 
Pgait) = P* ag {t), (31) 



-i(uo+ni 2 )t„-i(r+r 12 )t n / n \ , (£+jjj) -i(^ +o 12 )t 

,9lUj+ (rH40l 2 ' 



x e -i( r + r i2)t [2fi 12 e- r * sin(2fii 2 t) + T(l - e~ rt cos(20i 2 i))] p es (0) 
+ i 7 ^^re- < ^+ n ")*e-*( r+r »)* [2012(1 - e- r * cos(2tt 12 t)) 

- r e - rt sin(2r! 12 t)] p es (0), 
Pgs (t) = p* g (t). (32) 

The sixteen density matrix elements given by Eqs. (|123"|) to (j3"2)) completely solve the master equation (TiT)f describing 
the dynamics of the two-qubit system interacting with the vacuum bath. For the case of the Dicke model, where 
fij/Xo ~ * 0, i.e., when the interatomic separation is much smaller than the resonant wavelength, (dp/dt) a = and the 
solution of the effective three-level system can be extracted out of the above equations. The conditions under which 
the Dicke model is obtained is analogous to the case of collective decoherence for the case of two-qubit interaction 
with a bath via a quantum nondemolition interaction (QND) [25j . There it was found that for the case of interaction 
with a thermal (and also a vacuum) bath, the subspace spanned by {|ei,<? 2 ), |<?i,e 2 )} is a decoherence-free subspace, 
implying that the matrix elements Pe 1 .g 2 ;e 1 ,g 2 , Pgi,e 2 ;gi.e 2 Pe 1 ,g 2 -g 1 .e 2 and p gi ,e 2 -e 1 ,g 2 remain invariant inspite of the 
system's interaction with a bath. However, from Eqs. (|23| (|32l) it is clear that none of the matrix elements is 
invariant as a function of time, reflecting the greater complexity of the dissipative interaction. 

B. Squeezed thermal bath 

Here we consider the two-qubit dynamics resulting from an interaction with a squeezed thermal bath, i.e., make use 
of Eq. ([5]). The equations of the reduced density matrix ([5]) taken in the two-qubit dressed state basis (fT%j) are not 
all mutually coupled, but divide into four irreducible blocks A, B, C, D, thereby reducing the task from an evaluation 
of fifteen coupled linear differential equations to that of a maximum of four coupled equations. Thus we have: 



Block A : 

p ee (t) = -2T{N + l)p ee (t) + n{(T + T 12 ) Pss (t) + (T - r 12 )p aa (t)} + r 12 \M\ Pu (t), 

p ss (t) = -(T + T 12 ){-N+(l + 3N)p ss (t)- p ee (t)+Np aa (t) + \M\p u (t)}, 

p aa {t) = {r-T 12 ){N-(l + 3N)p aa (t)+p ee {t)-N Pss (t) + \M\p u (t)}, 

Pu (t) = 2|M|r 12 -4u;o|p S e|sin($ + x)-(2iV^^ (33) 

Here 

p u (t) = e i *p ge (t) + h.c, (34) 

Pge = \Pge\e lx and $ is as in Eq. (fT0|) . Also p gg (t) = 1 — p aa {t) — p ss (t) — Pee{t)- The Eqs. (|33|) give the dynamics of 
the population of the two-qubit system interacting with a squeezed thermal bath while the off-diagonal terms, given 
by the Blocks B, C and D are: 
Block B : 

Pes (t) = -%(uq - n 12 )) Pes (t) - i{(3r + r i2 ) + 2A>(2r + ri 2 )} Pes (t) + n{t + r 12 ) Psg (t) 

+ MYnPgsit) - M(T + T 12 )p se (t), 
p se (t) = p* es {t), 

Pgs (t) = -*(wo - n 12 )) Pgs (t) - -{(r + r 12 ) + 2A>(2r + r 12 )} Pgs (t) + (1 + n)(t + r 12 ) Pse (t) 

+ M*T 12 p es {t) - M*(T + T 12 )p sg (t), 
Ps B {t) = P* gs {t). (35) 

Block C : 

Pas (t) = i2tl 12 p as (t)-T(l + 2N)p as (t), 

Psa{t) = P * as {t). (36) 

Block D : 

Pea (t) = -»(wo + n 12 )) Pea {t) - -{(3r - r i2 ) + 2A>(2r - r 12 )} Pea {t) - n(t - v 12 ) Pag {t) 

+ MT 12 p ga (t)- M(T-T 12 )p ae {t), 

Pae(t) = P* ea (t), 

Pga (t) = i(tu Q - n 12 )) Pga (t) - -{(r - r i2 ) + 2A>(2r - r 12 )} Pga (t.) - (i + n)(t - r 12 ) Pae (t) 

+ M*T 12Pea (t) - M*(T - T 12 )p ag (t), 

Pag{t) = P* ga (t). (37) 

The Eqs. (|3"6"|) can be trivially solved to yield 

p as (t) = e^-^+^p as (0) 

p sa (t) = p* as (t), (38) 

while the blocks A, B and D consist of four linear coupled differential equations which can be written in matrix form 
as: 

P(t) = -QP{t) + W, (39) 

where Q is a time- independent 4x4 matrix and Pit), W are 4x1 column vectors. The Eq. ([39|) has the general 
solution 

P(t) = (Ve- Dt V^)P{0) + VD-\l - er m )V- x W. (40) 
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purity purity 




(a) (b) 



FIG. 1: Purity as a function of temperature T (in units where % = /cb = 1) for (a) the independent decoherence model, where 
Tij/Xo (|14|) is > 1 and (b) the collective decoherence model, where rtj/Xo (|14l) is w 0. Here with ri2 is the inter-qubit distance. 
The large-dashed, bold and dotted curves correspond to evolution time t = 1.0 and bath squeezing parameter (f9l llOfl r = —0.5, 
1.0 and 1.5, respectively. Here and in all the subsequent figures, the squeezing parameter $ (|10p is set equal to zero. Also cjo 
(1111) and the bath parameter F (|15l) . are set equal to 1.0 and 0.05, respectively. All the inter-qubit distances are defined on 
the scale of the resonant wavelength coming from the wavevector k ([7]l as a result of the position dependent couplings of the 
qubits with the bath. In figure (a) related to the independent decoherence model, kri2 is set equal to 1.5 while in figure (b) 
related to the collective decoherence model, fcri2 is set equal to 0.08. 



Here V is the vector composed of the eigenvectors of the matrix Q, while D is composed of its eigenvalues. We solve 
Eq. f4T)|) by numerically obtaining the eigenvalues and eigenvectors of the matrix Q for the Blocks A, B, C, and D. 

Figures Q](a)), (b) depict the behavior of purity, defined here as Tr(p 2 (t)) for p(t) as obtained in this subsection for 
the independent [k^.r^ ^ 0) and collective (fco-^ij - * 0) decoherence model, respectively, as a function of temperature 
T for an evolution time t and bath squeezing r ^\ [TO]). In all the figures in this article, we consider the initial state 
of one qubit in the excited state |ei) and the other in the ground state \g2), i-e., | ei } | <?2 ) and /x.fy (jTTJ) is equal to 
zero. It can be seen that with the increase in temperature, as also evolution time t and bath squeezing r, the system 
becomes more mixed and hence loses its purity. 



V. ENTANGLEMENT ANALYSIS 



In this section, we will study the development of entanglement in the two qubit system, both for the independent 
as well as the collective decoherence model interacting with a squeezed thermal bath. A well known measure of MSE 
is the concurrence [28[ defined as 

C =max(0, a/A-2- V^3 - a/aI), ( 41 ) 

where are the eigenvalues of the matrix 

R = pp, (42) 

with p = o y ® a y p*<jy ® a v and Oy is the usual Pauli matrix. C is zero for uncntanglcd states and one for maximally 
entangled states. Since the reduced dynamics of the two-qubit system was obtained in the dressed state basis (|18|) . 
which contains entangled states, in order to remove spurious entanglement coming from the basis, for the entanglement 
analysis we rotate the density operator back to a separable basis by means of a Hadamard transformation acting in 
the subspace spanned by {|a), |s)}, i.e., the tensor sum of a Hadamard in this subspace and an identity operation in 
the subspace spanned by {|e), \g}}, i.e., -ff( as ) © I{eg)- 

Here we study concurrence for the two-qubit system interacting with a squeezed vacuum bath. In figure (|2"|) 
concurrence is plotted for the initial state |ei) ^2) f° r both the independent as well as collective dynamics. Figure ([3]) 
depicts the behavior of concurrence for the same initial state with respect to the inter-qubit distance ri2- It is clearly 
seen that the buildup of entanglement is greater for the collective dynamics when compared to the independent one. 

Now we take up the issue of entanglement from the perspective of the PDF as in Eq. Q. In figures (0] (a)) 
and (b), we plot the weights u>i, 0J2, ^3 and ([3]) of the entanglement densities of the projection operators of the 
various subspaces which span the two qubit Hilbert space with respect to the evolution time t for the independent 
and collective decoherence models, respectively, for the case of an interaction with an unsqueezed vacuum bath. Since 
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FIG. 2: Concurrence C (|41[) as a function of time of evolution t. Figure (a) deals with the case of vacuum bath (T — r = 0), 
while figure (b) considers concurrence in the two-qubit system interacting with a squeezed thermal bath, for a temperature 
T — 1 and and bath squeezing parameter r ([9] I10[l equal to 0.1. In both the figures the bold curve depicts the collective 
decoherence model (kr\2 — 0.05), while the dashed curve represents the independent decoherence model (kr±2 = 1.1). In figure 
(b) for the given settings, the concurrence for the independent decoherence model is negligible and is thus not seen. 
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FIG. 3: Concurrence C (|41[) with respect to inter-qubit distance 7-12. Figure (a) deals with the case of vacuum bath (T = r = 0), 
while figure (b) considers concurrence in the two-qubit system interacting with a squeezed thermal bath, for T = 1, evolution 
time t = 1 and bath squeezing parameter r (fD TllOp equal to 0.1. In figure (a) the oscillatory behavior of concurrence is stronger 
in the collective decoherence regime, in comparison with the independent decoherence regime (fcri2 > 1). In figure (b), the 
effect of finite bath squeezing and T has the effect of diminishing the concurrence to a great extent in comparison to the vacuum 
bath case. Here the concurrence for the independent decoherence regime is negligible, in agreement with the previous figure. 



u)\ is the weight of the one dimensional projection, representing a pure state, and <x>4 that of the maximally mixed 
state with loi and L03 being intermediary, these plots depict the variation in the contribution of the various subspaces 
to the entanglement of the two-qubit system as t increases. From the figures it can be seen that in the case of the 
independent decoherence model, as depicted in figure (|3](a)), the weight 0J2 dominates the other weights and remains 
almost constant, while the remaining weights are much lower and their increase is very small compared to it. This 
is in contrast to the collective decoherence model, wherein we find the weight u>\ decreases while 102 increases with 
time. Since the weight lo\ is indicative of pure state entanglement, it is evident from the figures that the entanglement 
content in the collective decoherence model is higher than that in the independent decoherence model. Also since the 
weight o>2, representing the weight of the two-dimensional projection operator corresponding to the PDF Vi{£)i has 
a rich entanglement structure, we come to the conclusion that for the case of the two-qubit interaction with a vacuum 
bath, the entanglement is preserved in the system for a long time. In both the figures, the weight W4, indicative of a 
completely mixed state, is negligible over the time range considered. 

In figures ([5] (a)) and (b), we plot the weights lo\, L02, ^3 and UJ4 ^ of the entanglement densities of the projection 
operators of the various subspaces which span the two qubit Hilbert space with respect to T for the independent and 
collective decoherence models, respectively. These plots thus depict the variation in the contribution of the various 
subspaces to the entanglement of the two qubit system as T increases. From figure (0 (b)), we can see that in the 
collective decoherence regime, the weight Wi initially falls and then stabilizes around a finite value while U4 rises, but 
remains well below the value ofwi, indicating that for the collective decoherence model, under the given settings, the 
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FIG. 4: The weights ([3]) as a function of evolution time t, for the case of an interaction with an unsqueezed vacuum bath. Figure 
(a) refers to the independent decoherence model, with krvi = 1.5 and (b) the collective decoherence model, with kr\2 = 0.08. 
In both the figures, the bold curve corresponds to the weight uji, while the large-dashed, small-dashed and dotted curves 
correspond to the weights L02, UJ3 and u)4, respectively. In both the figures, the weight u>4 is negligible and hence is not seen. 

weights weights 




(a) (b) 



FIG. 5: The weights (0 as a function of temperature T, for the case of an interaction with a squeezed thermal bath for an 
evolution time t = 5 and bath squeezing parameter r I10|) equal to 0.5. Figure (a) refers to the independent decoherence 
model, with kr-12 = 1-5 and (b) the collective decoherence model, with kr-12 = 0.08. In both the figures, the bold curve 
corresponds to the weight ui, while the large-dashed, small-dashed and dotted curves correspond to the weights 0*2, u>3 and W4, 
respectively. 



system maintains a finite value of entanglement. This feature is not observed in the figure (JS](a)), where for the same 
settings, ui\ decreases to zero, while LO4 rises, thereby indicating a loss of purity and destruction of entanglement. 

As explained in Section II, the characterization of MSE for a two-qubit system involves the density function of four 
projection operators, IIi, II2, II3, II4, corresponding to one, two, three, and four dimensional projections, respectively. 
These will be represented here as V\(£), V2{£), V^{£) and Vi{£), respectively. As also discussed above, Va(£) would 
be universal for the two-qubit density matrices and would involve the Harr measure on SU(4) [38| . This is depicted 
in figure © and is common to all the two-qubit PDF of entanglement. 

Now we consider the Vi(£) and Vz{£) density functions for some representative states of the two qubit system, 
both for the independent as well as collective decoherence models. This enables us to compare the entanglement in 
the respective subspaces of the system Hilbert space. We also plot the full entanglement density function curve V(£ ) 
with respect to the entanglement £, at a particular time t. This will enable us to look at the contribution to the 
entanglement from the different projections. 

Figures ([7] (a)) and (b) depict the behavior of the density function Vii£) for the bath evolution time t = 5.0 
and T = for the independent and collective decoherence models, respectively. For these conditions, the value of 
concurrence ([4"Tj) is 0.17 for the case of the independent decoherence model and 0.42 for the collective model, depicting 
the greater entanglement content in the later compared to the former. This is also borne out by these figures. As shown 
in [3l|, the concurrence for a two dimensional projection is Cn 2 — (£max~ £c«sp)/2- From the figures |7|). it can be seen 
that the value of Cn 2 is greater for the case of figure (JTfb)) when compared to that of figure J7£a)). In all the figures 
related to the probability density functions of entanglement for the two-qubit system interacting with an unsqueezed 
vacuum bath, for the independent decoherence model, kr\2 is set equal to 1.5, while for the collective decoherence 
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model, kr\2 is set equal to 0.07. For the two-qubit interaction with a squeezed thermal bath, the independent 
decoherence model is parametrized as above, while for the collective decoherence model, kr 12 is set equal to 0.08. 
Figures (0(a)) and (b) exhibit the entanglement density function V^E) while the figures ([9] (a)) and (b) illustrate the 
behavior of the full density function V(£) for the independent and collective decoherence models, respectively and 
for the same parameters as above. The full entanglement density function is obtained by a weighted sum over all the 
contributions of the projection operators from the different subspaces ([3]). Here and in all the subsequent figures for 
the full PDF, the contribution from the one dimensional projection Hi, which is a delta function, is represented by a 
line of height equal to its weight ([3J and the point on the absicca is determined by its corresponding entanglement. 
The figures ([9j (a)) and (b) depict the rich entanglement structure present in the two-qubit system interacting with an 
unsqueezed vacuum bath. While the contribution in the figure ©(a)), for the independent decoherence model, comes 
primarily from the two dimesional projection, that for the collective decoherence model, as in figure (0(b)), comes 
from the one and two dimensional projections with the one dimensional projection being the dominant contributor. 
Since the one dimensional projection represents the pure state entanglement extant in the system, this clearly shows 
the greater entanglement content in the collective decoherence regime as compared to the independent one. The figures 
([TOl (a)) and (b) show the behavior of the full density function V(£) for the independent and collective decoherence 
models, respectively, for the bath evolution time t = 20.0 and T = 0. For these parameters, the value of concurrence 
(|4"Tj) is 0.32 for the case of the independent decoherence model and 0.54 for the collective model, depicting the greater 
entanglement content in the later compared to the former. By a comparison with the earlier case for an evolution 
time t = 5, is brought out the point that with time there is a build up of entanglement in the two-qubit system. 
The system is still possessed of a rich entanglement structure. For the independent decoherence model, the principal 
contribution still comes from the two dimensional contribution but in the collective decoherence regime, in contrast 
to the earlier case, the principal contribution has now shifted from the one to the two dimensional projection, thereby 
showing that the system gets more mixed with the passage of time. 

We now consider the nature of entanglement in the two-qubit system when it interacts with a squeezed thermal 
bath. Figures (|lll (a)) and (b) depict the behavior of the entanglement density function Vi(£) for the evolution time 
t = 1, T = 10 and bath squeezing parameter (HUTU} r = 0-5 for the independent and collective decoherence models, 
respectively. For these conditions, the value of concurrence (|4Tj) is 0, indicating a complete depletion of entanglement. 
This is partially borne out by the fact that for this case Cn 2 ~ 0. However, as seen from the figures (|T2l (a)) and 
(b) and also from the figures (fT3l (a)) and (b), the system still exhibits a rich entanglement structure. In figure (fT2l 
(a)), the parameter £j_ is « 0.48 while it is w 0.98 for the collective model as in (fT2l (b)). Since £± is evaluated from 
the state perpendicular to the third non-separable basis in the canonical basis of a three dimensional projection [31j . 
this clearly brings out the greater entanglement content in the collective model as compared to the independent case. 
In figure (TI31 (a)), the contribution to the full entanglement density for the independent model comes from the one, 
two and four dimensional projections with an approximately equal contribution from the one and four dimensional 
projections while in figure {131(b)), the contribution to the full entanglement density for the collective model comes 
from the one and four dimensional projections with the contribution from the one dimensional projection being the 
dominant one. Since the contribution from the one dimensional projection is an indicator of pure state entanglement 
in the system, this clearly brings out the entanglement content in the system, under the given conditions, even though 
concurrence is zero. Also the greater weight carried by the one dimensional projection in the collective model in 
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FIG. 7: The density function Vi(£) with respect to the entanglement £ for an evolution time t = 5.0, T = 0.0 and bath 
squeezing parameter r equal to 0. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence 
model. 



comparison with the independent one, brings out the greater entanglement content in the former. 

Figures (fill (a)) and (b) illustrate the full entanglement density for the evolution time t = 1, T = 10 and bath 
squeezing parameter Q [TU| r = 1.0 for the independent and collective decoherence models, respectively. For these 
conditions, the value of concurrence (|41j) is zero. However, the system still exhibits an entanglement structure, 
especially in the collective decoherence regime as seen in the figure (fill (b)) where the one, three and four dimensional 
projections contribute with the principal contribution coming from the one dimensional projection, thereby indicating 
the presence of pure state entanglement in the system. The corresponding case for the independent model, as in figure 
(j 141 (a)), has its principal contribution coming from the four dimensional projection representing a maximally mixed 
state. This clearly shows the richer entanglement content in the collective model in comparison to the independent 
one. Figures (fT51 (a)) and (b) illustrate the full entanglement density, for the same conditions as above but with 
zero bath squeezing, for the independent and collective decoherence models, respectively The entanglement in the 
system is seen to be greater for this case when compared to the previous cases with finite bath squeezing. The value 
of concurrence (|4T|) is zero for the independent model, while it is 0.17 for the collective one, bringing out the greater 
entanglement content in the later in comparison with the former. This is also borne out by the entanglement density 
function, where the contribution to the full density function from the one dimensional projection (indicative of pure 
state entanglement) is greater in the collective model when compared to the independent one. Also, in this case of 
interaction with an unsqueezed thermal bath, the full entanglement density function for the independent decoherence 
model exhibits greater entanglement when compared to the corresponding case of interaction with a squeezed thermal 
bath. This seems to indicate that for the two-qubit system, a finite bath squeezing is detrimental to the development 
of entanglement. Finally in figures (|16l (a)) and (b) is depicted the full entanglement density for the evolution time 
t = 5 , T = 10 and bath squeezing parameter ([3J [TUJ) r = 0.5 for the independent and collective decoherence models, 
respectively. As expected, with the increase in the exposure to the environment, indicated by the greater evolution 
time, the system becomes more mixed and hence loses entanglement. This is borne out by the fact that for these 
conditions, the value of concurrence (|41|) is zero. However, for the collective decoherence model, as shown in figure 
(|16l (b)). the full density function has contributions coming from the one and four dimensional projections, with the 
weight carried by the one dimensional projection greater than that by the four dimensional one bringing out the fact 
that pure state entanglement is still extant in the system under the given conditions. However, for the independent 
decoherence model, as in figure (|16l fa)), the principal contribution to the full entanglement density function comes 
from the four dimensional projection, indicating that the system is tending towards a maximally mixed state, thereby 
losing its entanglement. 

Thus we see that when compared to the interaction with a bath at T = 0, the finite T bath is detrimental to the 
generation of entanglement between the two-qubit system. An interesting connection of this can be made with the work 
presented in (39| , where the authors connected frustration with interaction strength and ground-state entanglement. 
They defined frustration as the non-commutativity of the local Hamiltonian (consisting of the single body terms) with 
the interaction Hamiltonian and argued that with the increase in interaction, ground-state entanglement in the system 
would increase. The case of dissipative S — R interaction with [Hs,Hsr] ^ would suggest a frustrated system in 
this nomenclature. For the interaction of the two-qubit system with a vacuum bath, entanglement is seen to rise with 
increase in the influence of the environment qualitatively agreeing with the above work. Finally entanglement would 
vanish because as a result of decoherence, the system would lose its quantum coherence and tend towards the classical 
regime. 
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FIG. 8: The density function P3(£) with respect to the entanglement £ for an evolution time t = 5.0, and zero T and bath 
squeezing. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence model. 
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FIG. 9: The full density function V{£ ) ([3]) with respect to the entanglement £ for an evolution time t = 5.0, and zero T and 
bath squeezing. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence model. 

VI. AN APPLICATION TO QUANTUM COMMUNICATION: QUANTUM REPEATERS 

The technique of entanglement purification [32j can be adapted for quantum communication over long distances, 
the key idea behind a quantum repeater [33| . The efficiency of quantum communication over long distances is reduced 
due to the effect of noise, which can be considered as a natural open system effect. For distances much longer than the 
coherence length of a noisy quantum channel, the fidelity of transmission is usually so low that standard purification 
methods are not applicable. In a quantum repeater set-up, the channel is divided into shorter segments that are 
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FIG. 10: The full density function V(£) (0) with respect to the entanglement £ for an evolution time t = 20.0, and zero T and 
bath squeezing. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence model. 
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FIG. 11: The density function Vi(£) with respect to the entanglement £ for an evolution time t — 1.0, T = 10.0 and bath 
squeezing parameter r equal to 0.5. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence 
model. 
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FIG. 12: The density function V^E) with respect to the entanglement £ for an evolution time t — 1.0, T = 10.0 and bath 
squeezing parameter r equal to 0.5. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence 
model. 



purified separately and then connected by the method of entanglement swapping. This method can be much more 
efficient than schemes based on quantum error correction, as it makes explicit use of two-way classical communication. 
The quantum repeater system allows entanglement purification over arbitrary long channels and tolerates errors on 
the percent level. It requires a polynomial overhead in time, and an overhead in local resources that grows only 
logarithmically with the length of the channel. 

We consider the effect of noise, introduced by imperfect local operations that constitute the protocols of entangle- 
ment swapping and purification , on such a compound channel, and how it can be kept below a certain threshold. 
The noise process studied is the one obtained from the two-qubit reduced dynamics via a dissipative system-reservoir 
interaction, studied above, instead of the depolarizing noise considered in [331 ] - Here we treat this problem in a sim- 
plified fashion, and study the applicability and efficiency of entanglement purification protocols in the situation of 
imperfect local operations. A treatment of partial teleportation of entanglement, conceptually equivalent to quantum 
entanglement swapping, in a noisy environment was made in [4fj| . 

A quantum repeater involves the two tasks of entanglement swapping, involving Bell-state measurements, and 
entanglement purification, involving CNOT gates. The Bell-state measurement may be equivalcntly replaced by a 
CNOT followed by a projective single-qubit measurement. In entanglement swapping, two distant parties initially not 
sharing entanglement with each other, but sharing entanglement separately with a third party, become entangled by 
virtue of a multi-partite measurement by the third party on the latter's two halves of entanglement. Entanglement 
purification involves two parties employing local operations and classical communication (LOCC) to improve the 
fidelity F of Einstein-Podolsky-Rosen (EPR) pairs they share, with respect to a maximally entangled state. The local 
operations involve two-qubit gates such as the CNOT operation, followed by single qubit measurement, and a possible 
discarding of an EPR pair. Provided F > 0.5, and at the cost of losing shared (impure) entanglement, the two parties 
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FIG. 13: The full density function V{£) with respect to the entanglement £ for an evolution time t = 1, T = 10.0 and 
bath squeezing parameter r equal to 0.5. Figure (a) refers to the independent decoherence model and (b) to the collective 
decoherence model. 




FIG. 14: The full density function V(£) © with respect to the entanglement £ for an evolution time t = 1, T = 10.0 and 
bath squeezing parameter r equal to 1.0. Figure (a) refers to the independent decoherence model and (b) to the collective 
decoherence model. 



can increase the fidelity of the remaining shared entanglement to 

pi = F 2 + [(1 - ^)/3] 2 (43) 

F 2 + [2F(l-F)/3] + (5/9)(l-F) 2 ' V ; 

where F and F' are, respectively, the input and output fidelities of the entanglement purification protocol proposed 
by Bennett et al. (32 ] . 

In the simplified scenario considered here, the output of the noisy CNOT is taken to be a mixed separable state, in 
place of a pure separable state that is obtained in the noiseless case. As a further simplification, in order to faciliate 
an easy connection with the purification protocol due to Bennett et al., this mixed state is assumed to be of the form: 

p(F) = F 2 \ U)U, I I + (1 - F 2 )| UXi, T I, (44) 

where | f ) (| j)) = \\) (| — |)). Thus, p(F) is a mixture in the two dimensional space spanned by {| |, J.), | J., f )}, 
parametrized by fidelity F, given by y/ (f, J T) ■!)■ The state is then the input to the purification protocol, 

whereby we obtain the output fidelity F' as a function of F. 

We depict in figures (fl"7l (a)) and (b) the modified 'purification loop', obtained by subjecting the noiseless loop to 
the above model of noise for the independent and collective decoherence regimes, respectively. In both figures we find 
that noise degrades the performance of the purification protocol. To evaluate its performance, it may be compared 
with the noiseless case, given by the dot-dashed curve in both figures. That this curve lies above the F = F' line in 
the closed range [0.5, 1] implies that fidelities above the minimum value F m i n = 0.5 can be corrected to the maximum 
value -Fmax = 1 by repeated application of purification. The degrading effect of noise can be seen in two ways: it 
introduces an off-set, whereby an input of F — 1 does not yield the same output; it restricts -Fmax to values less than 
1. 
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FIG. 15: The full density function V(£) with respect to the entanglement £ for an evolution time t = 1, T = 10.0 and 
bath squeezing parameter r equal to 0, i.e, we consider here an unsqueezed thermal bath. Figure (a) refers to the independent 
decoherence model and (b) to the collective decoherence model. 




FIG. 16: The full density function V{£) © with respect to the entanglement £ for an evolution time t = 5, T = 10.0 and 
bath squeezing parameter r equal to 0.5. Figure (a) refers to the independent decoherence model and (b) to the collective 
decoherence model. 



In figure (|17l (a)), comparison of the small and large dashed curves shows that increasing bath exposure time 
degrades the fidelity by suppressing F max (the point where a curve cuts the F = F' line from above) and increasing 
the off-set at F = 1. Although figure (|17l (b)) similarly shows the expected off-set due to noise, the surprising feature 
is that the curves for t = 3 and t = 20 show a slight lowering of output fidelity with increasing input fidelity. The 
matter is compounded by noting that at an intermediate time (t = 14), as depicted by the largc-dash-dot curve, the 
function F'(F) is monotonously increasing. This may be attributed to the fact that entanglement generated due to 
interaction with the bath shows a strong oscillatory behavior, as seen from figure © for the collective regime. 



VII. CONCLUSIONS 



Here we have analyzed the dynamics of entanglement in a two-qubit system interacting with its environment, taken 
to be in a general squeezed thermal state, via a dissipative S — R interaction. The analysis of the mixed state 
entanglement has been made using a measure involving a probability density function (PDF). 

The position dependent coupling of the qubits with the bath enabled a natural division of the dynamics into an 
independent and collective decoherence regime, where in the independent decoherence regime the qubits interact via 
localized S — R interactions, while in the collective regime the qubits are close enough to feel the bath collectively. 
The reduced dynamics revealed that in this case of dissipative S — R interaction, there is no decoherence-free subspace 
in contrast to that of a QND S — R interaction [25]. This can be understood in the following manner. As shown 
in [25|, the two-qubit QND reduced dynamics obeys the relation p a bif) = £ab(t)p a b(0), with the non-trivial aspect 
of the dynamics being that C represents, not a matrix, but a two-dimensional array, with the multiplication done 
element-wise. This gets translated into the spin-flip symmetry obeyed by the reduced dynamics. The emergence of 
a decoherence-free subspace could be attributed to the non-trivial symmetry obeyed by the system. An absence of a 
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FIG. 17: Purification loop for connecting and purifying EPR pairs. The noisy channel is modelled as a two-qubit dissipative 
interaction with an unsqueezed vacuum bath in the independent (kr a t > 1) (plot (a)) and collective (kr a b <C 1) (plot (b)) 
decoherence regime, with the input state given by p(F) in Eq. (|44|l . The bold line is the F — F' plot, the small and large 
dashed curves represent t — 3 and t — 10, respectively, while the dot-dashed curve is due to the noiseless Bennett et al. protocol. 
In (b), the larger dash-dot curve and the dotted curves represent, respectively, t = 14 and t = 20, respectively. 



decoherence-free subspace in the present case of dissipative interaction could be due to the lack of such a symmetry 
in the reduced dynamics. 

We studied the entanglement in the two-qubit system for different bath parameters, for both vacuum as well as 
a bath at finite T and squeezing and used the measure of mixed state entanglement involving a PDF as well as 
concurrence. It clearly emerged that entanglement generation is more efficient in the case of interaction with a bath 
at zero T than at finite T. It also appears that in this case of two-qubit dissipative S — R interaction, presence of finite 
bath squeezing does not help in entanglement generation. This is in contrast to expectations from the single-qubit 
dissipative interaction [36| where it was observed that for some bath parameters, the process of decoherence could be 
slowed down resulting in the preservation of quantum coherence for a longer time. It would be pertinent to point out 
here an interesting work (4l| where the authors studied in detail the evolution of entanglement between two oscillators 
coupled to the same environment. There the oscillator system was assumed to start from a general two- mode Gaussian 
state, which includes the two-mode squeezed state. Entanglement in the final state was considered as a resource whose 
origin could be the quantum resource available in the system, viz. squeezing in the initial oscillator system. It was 
shown that under certain regimes of the squeezing parameter, the environment could be used as a resource to extract 
entanglement, thereby highlighting the importance of squeezing in the study of entanglement generation. The work 
presented here differs from that in [4l|, apart from the obvious difference in the system (here a two-qubit system), 
in that the evolution considered here is Markovian while that in (4lj is that of quantum Brownian motion which is 
essentially non-Markovian. Also, here the squeezing originated from the bath while that in [f0 | was in the initial 
system state. It would be of interest to extend the work presented here to have a better understanding of the impact 
of squeezing, a natural resource, on the dynamic evolution of entanglement. 

The entanglement analysis via the PDF showed an advantage over concurrence, specially in the case of mixed state 
entanglement due to interaction with a finite T bath. While concurrence was zero for most of the cases considered, 
the PDF revealed the structure of entanglement still present in the two-qubit system bringing out the statistical and 
geometrical nature of the measure. Finally we made an application of the dissipative two-qubit reduced dynamics to 
a simplified model of a quantum repeater, which can be adapted for quantum communication over long distances. 
Thus this work along with [25[ makes a detailed study of the dynamics of two-qubit entanglement in the presence of 
generic, purely dephasing as well as dissipative, open quantum system effects. 
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